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Abstract 

Elk (Cervus canadensis ) historically are among the most widely distributed members of the 
deer family, occupying much of the United States, Canada, and northern Mexico. The natural 
distribution of this species decreased substantially in the early 20 th century, presumably resulting 
in the extirpation of populations in Texas. In the past 40 years, several herds of free-ranging 
elk have re-appeared in the Trans-Pecos region of Texas. For some herds, it is not known if 
the origin was: 1) the result of individuals that escaped from captive herds; 2) an expansion of 
previously transplanted individuals from South Dakota and Oregon into Texas; or 3) the result 
of natural emigrants from southeastern New Mexico into the Trans-Pecos region. The objective 
of this study was to use DNA sequences from the mitochondrial cytochrome-# gene and D-Loop 
region, in combination with nine microsatellite loci, to assess genetic divergence, relationships, 
and origin(s) of the contemporary elk herds in Texas. Findings of the mitochondrial sequence 
data depicted a high degree of relatedness among individuals throughout the sampling area; 
whereas, microsatellite data revealed differences in frequencies of alleles in the Glass Mountain 
populations of Texas compared to samples from South Dakota, New Mexico, and the Davis 
Mountains. Further, computer simulations of population genetic parameters based on the mi¬ 
crosatellite data supported a scenario depicting the origin of contemporary elk in Texas likely 
was the result of natural emigrants from New Mexico or descendants of previously introduced 
individuals from New Mexico. In addition, simulations did not detect evidence of a genetic 
bottleneck during the past 350 generations, indicating a long, shared history between Texas and 
New Mexico populations. 
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Introduction 


Prior to European settlement, wapiti or elk 
(Cervus canadensis ) were among the most widely 
distributed members of the deer family on the North 


American continent, numbering over 10 million indi¬ 
viduals (Davis 1940). The distribution historically was 
partitioned into eight recognized subspecies distributed 
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across northern portions of the North American grass¬ 
lands (modern day Canada and Alaska) southward to 
the southern edge of the Great Plains grasslands into 
modern day Mexico (Davis 1940; Hall 1981; O’Gara 
and Dundas 2002). Extensive sightings of elk were 
documented from 1600 to 1877 followed by drastic 
drop off in the population following in the 1880s (Davis 
1940). Due to anthropogenic events in the late 1800s, 
current distribution ranges and population sizes across 
all subspecies of elk were reduced to approximately 
10% of their historical North America geographic 
range and total population numbers were reduced to 
approximately 41,000 individuals in the late 1880s 
(Nowak 1999; Bryant and Maser 1982). By the 1900s, 
market hunting, sport hunting, habitat modification, and 
general agriculture practices further reduced overall 
population numbers and in some cases, subspecies such 
as C. c. merriami and C. c. canadensis were driven to 
extinction (Davis 1940; McCullough et al. 1996; To- 
weill and Thomas 2002). Bailey (1905) reported that 
elk had not been observed in the Trans-Pecos region 
for 20 years suggesting that elk, presumably, were 
extripated in Texas. 

Isolated populations currently are known from 
the Trans-Pecos and Panhandle regions of Texas. They 
include herds in the Davis, Glass, Guadalupe, Sierra 
Diablo, and Wylie Mountains, and in the northern Pan¬ 
handle (Witt 2008; Gill 2013; Schmidly and Bradley 
2016). The Trans-Pecos area elk population reached 
an estimated 400 individuals by the late 1930s and 
has gradually increased in recent years (Wright and 
Thompson 1935; Davis 1940; Genoways et al. 1979; 
Witt 2008; Coykendall 1990). In 2016, the population 
of elk in the Trans-Pecos region was estimated at ap¬ 
proximately 3,500 individuals (Schmidly and Bradley 
2016). The origin of these current populations has been 
highly debated for the past 30 years due to incomplete 
or uncertain records pertaining to historic introduction, 
importation for game ranches, and a possible influx 
of individuals into Texas from neighboring popula¬ 
tions in New Mexico—which themselves probably 
are products of introductions from Colorado (Ligon 
1927; New Mexico Department of Game and Fish 
1967; Findley 1975). 

At least three possible scenarios can explain the 
origin of contemporary populations of elk in the Trans- 
Pecos region of Texas. First, it is well-documented that 


individuals were introduced in the Trans-Pecos region 
several times by private parties (Davis 1940; Swepston 
1985). Further, many of the high-fence hunting facili¬ 
ties, located throughout the Trans-Pecos region, have 
introduced and bred captive-reared individuals to 
generate populations for hunting revenue. Several of 
these herds are routinely observed in the lower eleva¬ 
tions surrounding Alpine, Fort Davis, Fort Stockton, 
Marathon Basin, Stockton Plateau, and other foothill 
and grassland regions in the Trans-Pecos region. Under 
this scenario, few records are available and few follow 
up studies have been conducted to evaluate the success 
of said introductions. It is possible that at least some 
current populations of elk in the Trans-Pecos are the 
result of escapees that later founded local populations; 
although many of these cases are difficult to document. 

The second scenario involves situations where 
previous introductions were conducted to establish 
free-ranging populations instead of captive popula¬ 
tions. We have treated this scenario separate from 
introductions listed under the first scenario because of 
the timing (some nearly one hundred years ago), these 
introductions are accompanied by reliable sources of 
documentation and historical data, and in some cases, 
they are re-occurring studies to monitor the success of 
the introductions. One of the best-documented intro¬ 
ductions occurred in 1928, when Judge J. C. Hunter 
and his associates released 44 elk, collected from the 
Black Hills of South Dakota, into McKittrick Canyon 
located in the Guadalupe Mountains of Texas (Davis 
1940; Schmidly 2002). This herd steadily increased 
in number and reached a size of approximately 350 
individuals by the mid 1960s. Although this herd de¬ 
clined in the 1980s (McAlpine 1990), a viable popula¬ 
tion remains in this area. Additional reintroductions 
in the Davis, Eagle, Glass, Sierra Diablo, and Wylie 
Mountains occurred on several occasions with varying 
success and levels of documentation (Swepston 1985; 
Coykendall 1990; Gill 2013). For example, a large- 
scale introduction in 1988 involved the transplant of 48 
elk into the Davis Mountains and 51 elk into the Wylie 
Mountains from an eastern Oregon source population 
(Coykendall 1990). It is possible that descendants of 
these introductions in the Davis and Wylie Mountains 
spawned the populations now inhabiting the Trans- 
Pecos region (Eagle, Glass, and Guadalupe Mountains, 
as well as the intervening lowlands). 
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A third scenario that could explain the repopula¬ 
tion of elk in the Trans-Pecos region is from natural 
emigration from nearby New Mexico populations. The 
Trans-Pecos elk populations are connected to current 
populations of elk in New Mexico by several potential 
routes including the Sacramento, Delaware, and Gua¬ 
dalupe Mountains; or alternatively lowland crossings 
across foothills and plains regions. These natural cor¬ 
ridors could be used for immigration of elk into Texas 
from established herds occupying New Mexico. It is 
possible that populations throughout the current distri¬ 
bution are conforming to the resident species theory, 
which promotes the presence of population-specific 
allele frequencies, and would indicate that populations 
between the two states have been separate and are pur¬ 
suing their own evolutionary trajectories (Pohler et al. 
2014). However, if populations are freely traveling be¬ 
tween the border of the states and breeding, admixture 
would result in an inability to detect any type of genetic 
structure between the two states. It is worth noting that 
elk were extripated in New Mexico around 1908 and 


modem-day populations are a result of reintroduction 
efforts involving Rocky Mountain Elk (C. e. nelsoni ) 
from Colorado (Ligon 1927; New Mexico Department 
of Game and Fish 1967; Findley 1975). 

The goal of this study was to determine if any 
of these three possible scenarios (or combination of 
scenarios) reflects the most appropriate explanation for 
the elk located in the Trans-Pecos region. Therefore, 
DNA sequences from the mitochondrial cytochrome-# 
(Cyt- b) and D-loop control region, and nine microsatel¬ 
lite markers were examined due to their ability to detect 
differences at the species, subspecies, and population 
level (Meredith et al. 2005,2007; Fudt et al. 2004). To 
elucidate among these three scenarios, samples from 
the Davis and Glass mountain ranges were compared 
to samples collected throughout the range of C. c. 
nelsoni. Further, samples from South Dakota (C. c. 
canadensis) were included to serve as representatives 
of historic individuals that were transplanted into the 
Trans-Pecos region. 


Materials and Methods 


Sampling .—Samples were collected from free- 
ranging populations in the Glass (January 2008 and 
September/October 2009) and Davis (August 2014) 
Mountains (see Appendix). Ear clips and hair samples 
were obtained from 13 different individuals that were 
captured during a three-year ecological and move¬ 
ment study conducted on elk in this region (Pohler 
et al. 2014). Six samples from the Davis Mountains 
were acquired opportunistically from road kill and 
velvet sheds by Texas Parks and Wildlife Department 
(TPWD). The New Mexico Department of Game and 
Fish (NMDGF) and TPWD provided tissues (muscle, 
hair) and ear clip samples (2014/2015 hunting season) 
from neighboring populations potentially connected 
through possible emigration routes. Eight opportunistic 
samples (i.e., road kill/chronic wasting studies) were 
collected by NMDGF in 2013-2014 and added to the 
sampling scheme. South Dakota Game Fish and Parks 
provided four samples from the Black Hills in the fall 
of 2013. The Natural Science Research Faboratory 
(NSRF) Museum of Texas Tech University, provided 
two samples collected in southern New Mexico. In 
total, 38 samples (19 from Texas, 15 from New Mexico, 


and four from South Dakota) from 10 localities were 
examined in this study (Fig. 1). 

Specific data were recorded including locality and 
sex, and each individual was assigned a TK number, a 
unique identification number of the NSRF. No voucher 
specimens were collected due to the game species 
status, however, all possible tissues (muscle, hair, ear 
clips, and liver) were sampled and stored at -80° C. 

DNA sequencing .—Genomic DNA (nuclear and 
mitochondrial) was extracted following the protocols 
of the manufacture for a Qiagen DNeasy blood and 
tissue kit (Qiagen Inc., Valencia, California). Poly¬ 
merase chain reaction (PCR) was used to amplify the 
entire Cyt -b gene (1,140 base pairs). Amplification 
followed methods from Fudt et al. (2004) using prim¬ 
ers FGF766/765 (Bickham et al. 1995; Bickham et al. 
2004). D-loop primers (F0, Dl, El, and SO; Douzery 
and Randi, 1997) were used to amplify a 1,296 base pair 
region. PCR methods followed the standard HotStar 
Taq (Qiagen Inc., Valencia, California) amplification 
protocol. PCR amplifications for both mitochondrial 
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Figure 1. Map depicting the location of sampling sites. White circles denote 
collection localities of individuals examined in this study and numbers refer 
to the localities provided in the Appendix. 


regions were conducted in 25-pL reactions containing 
25 ng of DNA, 5 units of Taq, 100 pM of each primer, 
25 pM of MgCl„ 10 pM of Bovine serum albumin, 
10 pM of a nucleotide mixture and 5.8 pL of double 
distilled water (ddH 2 0). The thermal profile for both 
Cyt -b and D-loop was as follows: an initial denaturation 
of 94° C for 3 min, followed by 28 cycles of 94° C 
for 30 sec, 45° C for 30 sec, and 72° C for 40 sec, and 
a final extension at 72° C for 20 min. PCR products 
were purified with ExoSAP-IT (Affymetrix, Santa 
Clara, California). 

PCR Primers and ABI Prism Big Dye version 
3.1 terminator technology (Applied Biosystems, Fos¬ 
ter City, California) were used to cycle sequence the 
amplified products. Cycle sequencing reactions were 


purified using Sephadex G50 (Amersham Pharmacia 
Biotech, Piscataway, New Jersey) cleanup protocols 
and were analyzed with an ABI 3100-Avant automated 
sequencer. Resulting sequences were aligned and 
proofed using Sequencher 4.0 or 4.1.2 software (Gene 
Codes, Ann Arbor, Michigan); chromatograms were 
examined to verily all base calls. All Cytb sequences 
obtained in this study were deposited in GenBank and 
are listed in the Appendix. 

Microsatellite amplification .—Nine microsatel¬ 
lites were selected from previous studies based on 
their effectiveness in detecting population-level allelic 
frequency differences in elk. Primers t-28, t-108, t-158, 
t-193, t-501, t-107, t-273, t-115, and t-507 (Meredith et 
al. 2005,2007; Table 1) were amplified in 33 individu- 
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Table 1. The primers and dye tags used in this study are indicated below. Standard allelic and genotypic informa¬ 
tion is included as follows: number of alleles, size range, observed heterozygosity (H o ), expected heterozygosity 
(H e ), and polymorphic informative content (PIC). 


Primer 

5’ dye Label 

No. Alleles 

Size Range 
(base pairs) 

H 

O 

H 

e 

PIC 

T26 

6FAM 

10 

328-398 

0.625 

0.802 

0.765 

T107 

6FAM 

15 

235-259 

0.441 

0.719 

0.695 

T108 

6FAM 

11 

136-181 

0.611 

0.879 

0.854 

T115 

6FAM 

12 

180-192 

0.742 

0.807 

0.775 

T156 

VIC 

22 

143-249 

0.813 

0.932 

0.911 

T193 

6FAM 

5 

184-220 

0.733 

0.727 

0.667 

T501 

NED 

15 

252-290 

0.531 

0.767 

0.729 

T507 

NED 

13 

148-202 

0.647 

0.862 

0.837 


als using PCR methods similar to Meredith (2007). 
The thermal profile includes a heating phase at 80° 
C, followed by 95° C for 10 minutes, then 40 cycles 
each of 45 second denaturation at 95° C, a 90 second 
annealing step from 48° C to 64° C (based on melting 
point for each primer) and 60 second elongation at 72° 
C with a final extension for 10 minutes at 72° C. Allele 
separation included 0.5 pi of 400 or 500HD RX size 
standard (Applied Biosystems Inc., Waltham, MA), 
8.5 pi of Hi-Di Formamide (Applied Biosystems Inc., 
Waltham, MA), and 1 pi of PCR product. GeneMapper 
software (version 4.0; Applied Biosystems Inc., Foster 
City, California) was used to determine fragment size, 
and ultimately allelic composition. Alleles were binned 
(averaged) to account for marginal error in determining 
fragment size. 

Data analyses. —Gene and nucleotide diversity 
were estimated using MEGA v6 (Tamura et al. 2013) 
for all individuals included in the study. The program 
MUSCLE (Edgar 2004) within the MEGA v6 soft¬ 
ware package (Tamura et al. 2013) was used to align 
30 sequences of both Cyt- b gene and D-loop region 
to outgroup and reference sequences of C. elaphus , 
Odocoileus virginianus, and Odocoileus hemionus. 

All full-length Cyt- b sequences were uploaded 
into Popart (Leigh and Bryant 2015). A parsimony 
network was constructed to determine the number of 


haplotypes throughout the distribution (Clement et al. 
2000). Results were then mapped over the distribu¬ 
tion/sample locality map of Texas, New Mexico, and 
South Dakota. 

To infer the evolutionary relationships of individ¬ 
uals, Bayesian phylogenetic analyses were conducted 
on the Cyt- b, D-loop, and combined dataset of both 
mitochondrial genes in MrBayes 3.1.2 (Ronquist and 
Huelsenbeck 2003). Cyt- b gene and D-loop sequences 
were aligned and partitioned by gene and Cyt- b se¬ 
quences were further partitioned by codon position. 
Model selection using MrModeltest 2.3 (Naylander 
et al. 2004) selected GTR +1 + Y model of nucleotide 
substitution for both genes; values for model parameters 
were not defined a priori but were treated as unknown 
variables with uniform priors. Bayesian analyses were 
performed using 4 simultaneous Markov chains at a 
“temperature” of 0.02, with random, unconstrained 
starting trees. Each chain ran for 20,000,000 gen¬ 
erations with trees sampled every 1,000 generations. 
Three independent MrBayes analyses were conducted 
to ensure that each run produced similar stationary 
likelihood values. The first 1,000,000 trees of each 
run were discarded as burn-in and the remaining trees 
were used to calculate the posterior probabilities and 
50% majority rule consensus tree. Clade probabilities 
of > 0.95 were used to indicate nodal support. 
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Microsatellite analysis .— Structure version 2.3.4 
(Pritchard et al. 2000) was used to estimate genetic 
population structure within the three sampling areas 
using the multi-locus genotype data. Initially, K (hy¬ 
pothesized population number) was set to 3 to avoid 
underestimation of the number populations represented 
in the study; assuming that K would equal either: all in¬ 
dividuals belong to a single population (K = 1), that the 
individuals belong to two populations (K = 2), and that 
the individuals belong to three or more populations (K 
= 3). These tests were then replicated two times using 
Admixture and no Admixture. Structure parameters: 
bum-in length = 1,000,000 Monte Carlo Markov chain 
repetitions = 10,000,000, K = 3, with ten iterations at 
each K value. The allelic dataset was uploaded into 
Structure, analyzed and the output data were compiled 
and compressed to be loaded into Structure Harvester 
to determine the most appropriate value of K (Earl et 
al. 2012). Highest value obtained for AK was used to 
determine the estimated populations of both runs, and 
the appropriate file of population assignment. 

The program GenAlEx version 6.501 (Peakall 
and Smouse 2006) was used to estimate allele frequen¬ 
cies as well as observed and expected heterozygosity 
(H o and H e ). The theoretical populations defined by 
the harvester analysis were used as a priori groupings 
for Analysis of Molecular Variance (AMOVA), and 
Principal Coordinate Analysis (PCoA). Probability 
of Identity analyses (probability of encountering 2 
individuals with identical genotypes) was calculated 
within GenAlEx. Polymorphic information content 
(PIC) was estimated using CERVUS version 2.0 (Mar¬ 
shall et al. 1998). 

Six DIY Approximate Bayesian Computation 
simulations (DIYABC, Cornuet et al. 2014) were de¬ 


signed to examine the three scenarios used to explain 
the possible origin of contemporary elk in Texas. Two 
simulations per scenario (one inferring a population 
bottleneck and one without a bottleneck) were run 
to more thoroughly assess potential variables associ¬ 
ated with the three possible scenarios. The DIYABC 
simulations were as follows: 1) fragmentation of Texas 
populations from the historic panmictic distribution 
(sensu early 1800 North America population) with a 
genetic bottleneck; 2) successful establishment of South 
Dakota elk into the Guadalupe Mountains of Texas 
(in approximately 1928) with a genetic bottleneck; 3) 
natural immigration of elk into Texas (approximately 
1930’s) from naturally occurring herds in New Mexico 
with a genetic bottleneck; 4) fragmentation of Texas 
populations from the historic panmictic distribution 
without a bottleneck {sensu 1800 North America 
population); 5) successful dispersal of elk into Texas 
from the South Dakota individuals that were introduced 
into the Guadalupe Mountains without a bottleneck 
(approximately 1928); and 6) natural immigration of 
elk into Texas (approximately 1930’s) from established 
herds in New Mexico without a bottleneck. 

Parameters of the DIYABC simulations (n = 
population size; t = time in generations) were based on 
historical data from Gill (2013) with a generation time 
of two years. The effective population size used in the 
simulation (n = 378) was estimated from the LDNE 
program (Waples and Do 2008). For simulations, 
parameters incorporated included: current popula¬ 
tion size, population size during inferred bottleneck, 
duration of bottleneck, pre-bottleneck population size, 
historic population size, and estimates of generations. 
Bayesian posterior probabilities were then used to 
compare the simulated and observed data. 


Results 


Thirty Cyt -b sequences were examined and 
possessed the following nucleotide composition: T = 
29.1%, C = 27%, A = 30.8%, and G = 13.1%. Thirty- 
one sites were determined to be variable and six sites 
were phylogenetically informative. The D-loop dataset 
consisted of 33 sequences and the nucleotide sequence 
composition as follows: T = 30.3%, C = 22.5%, A = 
31.5%, and G = 15.8%. Eighty-four sites were deter¬ 


mined as being variable, 74 of these sites were the result 
of insertions or deletions, and eight where identified as 
standard phylogenetically informative sites. 

Cyt -b and D-loop sequences were analyzed 
individually and then concatenated into a single da¬ 
taset (n = 30). All three Bayesian analyses produced 
similar topologies: consequently, only the results of 
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the concatenated analysis are discussed herein (Fig. 2). 
Nodal support was recovered for the separation of the 
European red deer (C. elaphus) and elk (C. canaden¬ 
sis ); however all samples of elk were contained with 
a single clade (albeit no differentiation of populations 
was supported). Genetic divergence values (Kimura-2 
parameter; Kimura 1980) indicated that individuals that 
comprised clades I and II differed by 0.38% and all 
individuals differed from European red deer by 5.73%. 

The microsatellite dataset (33 individuals) pos¬ 
sessed allelic richness values ranging from 5 (locus 


T-195) to 22 alleles (locus T-156) with an average 
of 12.88 (Table 1). Observed heterozygosity values 
(Table 1) ranged from 0.441 (locus T-107) to 0.813 
(locus T-156). Probability of Identity was 4.0 X 10' 8 
(1 in 25,000,000 individuals). In general, the results of 
the Structure Harvester analysis indicated that samples 
could be separated into two populations (Fig. 3). Struc¬ 
ture results indicate that Population 1 contained 13 
individuals from the Glass Mountains (n = 11), Davis 
Mountains (n = 1), and New Mexico (n = 1) with a 
posterior probability > 0.50; and that Population 2 was 
comprised of 14 individuals from the Glass Mountains 


■ Odocoileus he mi onus 
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Figure 2. Bayesian tree generated from the concatenated sequence dataset (Cyt -b and D-loop sequences). Roman 
numerals indicate clades as referred to in the text. Asterisks (*) indicate a posterior clade probability of 95% or 
greater. Individual labels are Glass Mountains (GM), Davis Mountains (D), New Mexico (NM), and South Dakota 
(SD). Locality and specimen data are provided in the Appendix. 
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(n =2), New Mexico (n = 8), and South Dakota (n = 4) 
with a posterior probability > 0.50. 

Structure population designations were used 
for the subsequent AMOVA and PCoA analyses. The 
AMOVA partitioned the genetic variation into 4% 
among populations and 96% within populations. The 
overall fixation index (F st population differentiation 
due to genetic structure) between Populations 1 and 2 
was 0.046 implying a panmictic population. Principal 
coordinate analysis (PCoA) was used to visualize re¬ 
lationships of inferred populations across the first two 
axes of variation (PCA 1 - 21.3% and PC A 2 - 12.2%; 
Fig. 4). Individuals from inferred Populations 1 and 
2 generally segregated primarily based on the second 
axis of variation. 

Microsatellite data processed in DIYABC showed 
support for one of the six simulations (natural immi¬ 
gration of elk into Texas from naturalized and native 
herds in New Mexico - without a bottleneck) was sup¬ 
ported (posterior probability value = 0.99) (Table 2). 
All other simulations resulted in posterior probability 
values < 0.45 leading to the dismissal of those simula¬ 
tion conditions as possible explanations for the origin 
of contemporary elk in Texas. Further, simulations did 
not detect evidence of a genetic bottleneck during the 
past 350 generations indicating a long, shared history 
between Texas and New Mexico populations. 
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Figure 4. A) Assignment of populations as determined from the principal coordinate analysis. 
Coordinates 1 and 2 explain 21.3% and 12.2% of the variation of the dataset, respectively. 
Groups are assigned based on designation by the Structure analysis. Individual labels are Glass 
Mountains (GM), Davis Mountains (D), New Mexico (NM), and South Dakota (SD). Black 
diamonds represent Population 1 and grey circles indicate samples assigned to Population 2 by 
Structure Harvester. B) Locus effect on the principle component analysis (PC 1 to PC2) depicted 
in Figure 4A. Direction and distance away from the origin (0.000) indicates the magnitude of 
variation contributed by each locus effecting the group association in Panel A. 
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Discussion 


Analyses of the concatenated dataset contain¬ 
ing the Cyt -b and D-loop sequences did not separate 
samples based on a recognizable geographic pattern 
(Figure 2); instead individuals from all major sampling 
regions (South Dakota, New Mexico, Davis Moun¬ 
tains, and Glass Mountains) were contained within 
both clades I and II. Strong support was obtained at 
the base of clades I and II (posterior probability value 
of 1.00); however, the genetic divergence between the 
two clades was small (0.35%). The lack of topological 
resolution from mitochondrial sequences is similar to 
that reported in other genetic studies of intraspecific 
relationships among elk (Ludt et al. 2004). It appears 
that these genetic markers are not evolving at a suf¬ 
ficient rate to distinguish among populations. 

Analysis of microsatellite data detected allele 
frequency differences among individuals sampled 
from the Trans-Pecos populations. Allelic richness 
values ranged from 5 (locus T-195) to 22 alleles (locus 
T-156) with an average of 12.88 (Table 1). Meredith 
et al. (2007) reported reduced microsatellite variation 
in Tule elk given the severe bottleneck effects in the 
late 1800s. Many of the microsatellite markers re¬ 
ported in their study possessed fewer than five alleles 
(lowest reported in this study was five with average of 
7.3). Low levels of genetic variability in the Tule elk 
was suggested to have occurred as a function of low 
numbers of founders rather than insufficient sampling 
(average number of alleles = 3.2). The genetic unique¬ 
ness of the Tule elk (Meredith et al. 2007) resulted 
from lack of genetic variation, not from novel genetic 
variability when compared against Rocky mountain elk 
(C. c. nelsoni ) and Roosevelt elk (C. c. roosevelti). In 
contrast, the current study detected a different pattern 
in allelic richness (20 unique alleles in Population 1) 
throughout the distribution of C. c. nelsoni with high 
levels of genetic diversity. 

Results of the Structure Harvester analysis indi¬ 
cated that the most appropriate partitioning of genetic 
differences is best explained by recognizing two popu¬ 
lations. The first population contained a total of 13 
individuals from the Glass Mountains (n = 11), Davis 
Mountains (n = 1), and New Mexico (n = 1), whereas 
the second population contained 14 individuals from the 


Glass Mountains (n =2), New Mexico (n = 8), and South 
Dakota (n = 4). When comparing the Glass Mountains 
(Locality 10) to the Davis Mountains (Locality 9), 
the two populations shared 73% of alleles with 27% 
unique to the Glass Mountains (20 alleles). The Davis 
Mountains individuals (Locality 9) were then compared 
to the individuals sampled closest in New Mexico (Lo¬ 
cality 8), which shared 88% of their alleles with 12% 
being unique to Locality 8. The results indicated the 
contemporary Texas herds could have originated from 
a New Mexico population dispersing into the Davis, 
Guadalupe, and Glass Mountains. Alternatively, the 
introduction of elk into the Glass Mountains in 1944 
(one bull and two cows, Swepston 1985) followed by 
gene flow from this “founder population” to populations 
in the Davis Mountains and New Mexico populations 
may explain the shared genotypes of these popula¬ 
tions and the presence of unique alleles in the Glass 
Mountains. This indicates that either elk in the Davis 
Mountains genetically are closer to elk in New Mexico, 
or alternatively there recently has been recent gene flow 
between the two populations. Further, the individuals 
from South Dakota did not form an independent group 
from either the Texas or New Mexico samples. 

The combined genetic data {Cyt- b, D-loop, and 
microsatellites) and computer simulations were used 
to evaluate the three major scenarios concerning the 
reappearance of elk in Texas. Both population Structure 
analysis and simulations indicated that the scenarios 
involving reintroduced or escapees from high fence 
ranches (Scenario 1), and the reintroduction of elk into 
Texas from South Dakota (Scenario 2) were not sup¬ 
ported. However, these results supported simulation 
6 (without an inferred bottleneck, posterior probability 
= 0.99), which modeled the immigration of elk into 
Texas from New Mexico without a bottleneck. The fact 
that simulations did not detect evidence of a genetic 
bottleneck during the past 350 generations depicts a 
long, shared history between Texas and New Mexico 
populations. 

Differences in the allelic abundance in elk 
populations in Texas support the recent and natural 
immigration of individuals from New Mexico or the 
survival of genotypes from past introductions involv- 
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ing individuals from New Mexico populations. It is 
possible that elk emigrated from either the Guadalupe 
Mountains or surrounding foothills into the northern 
regions of the Trans-Pecos via the Davis, Sacramento, 
and Delaware Mountain ranges. If descendants of 
the introduced individuals transplanted into the Davis 
Mountain region (from eastern Oregon) contributed 
significantly to the gene pool, then the contempo¬ 
rary samples from the Davis Mountain should not 
have shown a genetic affinity to the Glass Mountain 
population or to the individuals sampled from New 
Mexico; instead they should genetically be different. 
Unfortunately, we were not able to included samples 


from eastern Oregon into this study; consequently 
that scenario cannot be rejected until such samples are 
examined. Consequently, the scenario that Texas elk 
are a natural population and a product of emigration 
from herds in southern New Mexico and subsequent 
establishment of viable populations remains a viable 
scenario. More thorough sampling of elk populations 
throughout the Trans-Pecos area and southeastern New 
Mexico is needed to determine genetic boundaries and 
possible population subdivsion. Further, samples from 
Colorado and Oregon are needed to establish genotypes 
for comparison to Texas populations. 
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Appendix 

Specimens examined in the DNA sequencing and microsateliite portions of this study. For each specimen, the 
collection locality, specimen identification number (SD = South Dakota, NM = New Mexico, D = Davis Moun¬ 
tains, and GM = Glass Mountains), museum identification number (TK), and GenBank accession number are 
provided in parentheses. GenBank accession numbers for cytochrome-6 sequences are to the left of the slash 
and D-loop accession numbers are to the right. For three samples (TK254263, TK254264, and TK174564) only 
microsatellite data were obtained. Localities refer to those shown on Figure 1. 

Cervus canadensis .—Locality 1. South Dakota: Pennington County; Custer State Park (SD1: TK254268, 
KY751417/KY751473 and SD2: TK254269, KY751416/KY751474), Deerfield (SD3: TK254266, KY751419/ 
KY751471), and Jewel Cave (SD4: TK254267, KY751418/KY751472). 

Locality 2. New Mexico: Coflax County; 13 mi SE of Raton (NM9: TK254258, KY751420/KY751470). 

Locality 3. New Mexico: Rio Arriba County; 3 miles south of Valle Grando Peak (NM4:TK254254, KY751425/ 
KY751466). 

Locality 4. New Mexico: McKinley County; Marquez Wildlife Management Area (NM2: TK254252, KY751427/ 
KY751464). 

Locality 5. New Mexico: McKinley County (NM3: TK254253, KY751426/KY751465). 

Locality 6. New Mexico: Catron County; Pelona Mountain (NM5: TK254255, KY751424/KY751467; NM6: 
TK254256, KY751423/KY751468; andNM7: TK254257, KY751422/KY751469). 

Locality 7. New Mexico: Otero County; 8 miles W. Sacramento (NM8: TK123174, KY751421/KY751450). 

Locality 8. New Mexico: Chaves County; Benito canyon in the White Mountain Wilderness (NM1: TK254251, 
KY751428/KY751463). 

Locality 9. Texas: Jeff Davis County; Davis Mountains (Dl: TK254259, KY751433/KY751459; D2: TK254260, 
KY751432/KY751460; D3: TK254261, KY751431/KY751461; D4: TK254262, KY751430/KY751462; D5: 
TK254263, NA/NA; D6: TK254264, NA/NA; and D7: TK254265, KY751429/KY751475). 

Locality 10. Texas: Brewster County; Glass Mountains (GM1: TK254270, KY751438/KY751451; GM2: 
TK254271, KY751437/KY751455; GM3: TK254272, KY751435/KY751456; GM4: TK254273, KY751436/ 
KY751458; GM5: TK254274, KY751439/KY751449; GM6: TK254275, KY751440/KY751454; GM7: 
TK254276, KY751434/KY751448; GM8: TK174564, NA/NA; GM9: TK254277, KY751443/KY751452; GM10: 
TK254278, KY751442/KY751447; GM11: TK254279, KY751444/KY751446; GM12: TK254280, KY751441/ 
KY751453; and GM13: TK254281, KY751445/KY751476). 

Cervus elaphus .—AB021099 

Odocoileus hemionus .—HM222707 


Odocoileus virginianus. —DQ379370 
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